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Abstract 

Using the nonequilibrium Green's function method, we theoretically study the Andreev reflection(AR) in a 
four-terminal Aharonov-Bohm interferometer containing a coupled double quantum dot with the Rashba spin-orbit 
interaction (RSOl) and the coherent indirect coupling via two ferromagnetic leads. When two ferromagnetic 
electrodes are in the parallel configuration, the spin-up conductance is equal to the spin-down conductance due to 
the absence of the RSOl. However, forthe antiparallel alignment, the spin-polarized AR occurs resulting from the 
crossed AR (CAR) and the RSOl. The effects of the coherent indirect coupling, RSOl, and magnetic flux on the 
Andreev-reflected tunneling magnetoresistance are analyzed at length. The spin-related current is calculated, and a 
distinct swap effect emerges. Furthermore, the pure spin current can be generated due to the CAR when two 
ferromagnets become two half metals. It is found that the strong RSOl and the large indirect coupling are in favor of 
the CAR and the production of the strong spin current. The properties of the spin-related current are tunable in terms 
of the external parameters. Our results offer new ways to manipulate the spin-dependent transport. 
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Background 

A quantum dot (QD) is an artificially low-dimensional 
structure that can be filled with electrons (or holes). Two 
or more QDs can be coupled to form multiple-QD sys- 
tems (i.e., artificial molecules). Because the degrees of 
freedom of the QDs are well controllable, it is possi- 
ble to add or remove the electrons in the QDs, and the 
QD system can be coupled via tunnel barriers to elec- 
trodes, in which electrons can be exchanged. Accordingly, 
the artificial molecule provides an excellent model sys- 
tem in which the thorough investigation of quantum 
many-body properties in a confined geometry can be per- 
formed [1-6]. Among the various multiple-QD systems, 
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an Aharonov-Bohm (AB) interferometer containing dou- 
ble QDs (DQDs) is of particular interest and importance, 
in which two QDs are embedded in the opposite arms of 
the AB ring, respectively, and they are coupled to each 
other via barrier tunneling. As a tunable two-level sys- 
tem, the parallel DQD system that can become one of 
the promising candidates for the quantum bit in quantum 
computation has received more attention [7-20]. How- 
ever, in an actual DQD system, the coherent indirect 
coupling between two QDs via a reservoir is very essen- 
tial. Kubo et al. introduced the parameter a characterizing 
the indirect coupling strength, and Gurvitz also indicated 
the fundamentality of the sign of the coherent indirect 
coupling parameter [21,22]. Kubo et al investigated the 
pseudospin Kondo effect in a lateral DQD system using 
the slave-boson mean-field method and found that the 
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exotic pseudospin Kondo effect occurs when a coherent 
indirect coupling is presented through the common reser- 
voirs [23]. Recently, Kubo and co-workers calculated the 
shot noise and Kondo effect in a DQD structure with the 
coherent indirect coupling. Their results demonstrate that 
the coherent indirect coupling can generate a novel anti- 
ferromagnetic exchange phenomenon [24]. Trocha and 
Barnas studied theoretically the spin-dependent transport 
through a DQD coupled to ferromagnetic leads. They 
observed that the Fano antiresonance of the linear con- 
ductance relies on the sign of the indirect coupling in 
the nondiagonal coupling elements [8]. Furthermore, the 
transport properties of a DQD system has been con- 
sidered in the orbital Kondo regime. That the Kondo 
temperature and Kondo resonances are susceptible to the 
coherent indirect coupling parameter is also revealed [25]. 
In addition, if a QD is formed in a semiconductor two- 
dimensional electron gas structure without the inversion 
symmetry in the growth direction, the Rashba spin-orbit 
interaction (RSOI) will emerge, and the RSOI can induce 
the spin-related phase factor in the tunneling matrix ele- 
ments and the spin-flip effect. The RSOI results from a 
relativistic effect at the low speed limit, and it can cou- 
ple the electron spin to its orbital motion, thus providing 
a possible way to control the spin degree of freedom by 
means of an external electric field. As a consequence, the 
coherent indirect coupling and the RSOI make the quan- 
tum transport through the QD systems rich and varied 
[26-30]. 

On the other side, the subgap transport through het- 
erostructures with nano-objects (such as QDs, molecules, 
nanowires, etc.) coupled to one conductor and another 
superconducting lead has attracted a great deal of atten- 
tion over the past years due to the fundamental physics 
and its potential applications [31-35]. Andreev reflection 
(AR) usually occurs in the hybrid systems, in which two 
electrons with opposite spins enter the superconductor 
from the normal metal region, leading to the formation 
of a Cooper pair in the superconducting region [36-38]. 
In comparison with the standard mechanism of nor- 
mal AR, the crossed AR (CAR) is a nonlocal dynamics 
process which occurs at the contact between a super- 
conductor and two normal leads, where two subgap 
electrons from different metals enter into the supercon- 
ductor and generate a Cooper pair there [39-42]. AR 
(or CAR) in nanoscopic heterostructures gives rise to a 
rich subgap structure in the current-voltage character- 
istics. Accordingly, understanding the AR and CAR has 
attracted theoretical and experimental attention mainly 
because the AR (or CAR) may create the entangled elec- 
trons in a solid-state device, and CAR can be readily 
probed by spin selection using ferromagnetic electrodes. 
This approach is almost unrealized for entangler devices, 
since projecting the spin will cause the destruction of 



entanglement [43]. Based on the CAR, the controlled 
Cooper pair splitting has been realized in terms of a 
two-quantum dot Y-junction [44], which opens a possi- 
ble route towards a test of the Einstein-Podolsky-Rosen 
(EPR) paradox and Bell inequalities in solid-state sys- 
tems. Herrmann et al. used carbon nanotube DQD 
as Cooper pair beam splitters and realized the quan- 
tum optic-like experiments with spin-entangled electrons 
[45]. These results show that the CAR has an impor- 
tant application in testing a fundamental property of 
quantum mechanics. 

To our knowledge, the AR in the DQD with a maxi- 
mum coupling I a I = 1 has been studied widely. How- 
ever, the quantum transport through a four-terminal AB 
interferometer including a DQD in the presence of the 
AR, the coherent indirect coupling, and RSOI is less 
explored. Motivated by recent theoretical and experimen- 
tal advances in the DQD systems [7-10,13,15,16,19,21- 
25,44,45], one may expect that the interplay of the coher- 
ent indirect coupling and the RSOI in the presence of 
the AR will add new physics to hybrid quantum sys- 
tems, which may have practical applications for future 
spintronics. Consequently, we investigate the AR in the 
above-mentioned system in this paper. It is found that 
the RSOI and a nonzero coherent indirect coupling cause 
the spin-polarized AR when the polarizations of two fer- 
romagnetic leads are parallel, but for antiparallel (AP) 
arrangement of the polarizations of two ferromagnetic 
leads, the CAR can contribute the spin-polarized AR con- 
ductance. We note that the convex shape of the Andreev- 
reflected tunneling magnetoresistance (ARTMR) versus 
the magnetic flux relies on the sign of the coherent indi- 
rect coupling parameter, and there are extreme values 
in the plot of the ARTMR versus the coherent indi- 
rect coupling parameter. Even the negative ARTMR also 
occurs. This is a spin valve effect in the AR process. 
It is interesting to note that the sign of the coherent 
indirect coupling parameter leads to the swap effect in 
the spin-polarized current plot, and the pure spin cur- 
rent can be produced when two ferromagnetic leads are 
fully polarized. The spin-dependent AR current can be 
controlled by means of the gate voltage, RSOI, mag- 
netic flux, and so on. These results provide the ways to 
manipulate the spin-dependent transport by means of the 
system parameters. 

Methods 

We consider a hybrid four-terminal AB interferome- 
ter including a parallel DQD coupled to two ferro- 
magnetic reservoirs and two superconductors as shown 
in Figure 1. The system is described by the following 
Hamiltonian: 

H = Hy +Hs+ Hdqd + Hj, (1) 
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Figure 1 Schematic diagram of a four-terminal AB 
Interferometer (color on line). The AB interferometer contains a 
coupled DQD with magnetic flux applied perpendicular to rings. 



superconductors is also assumed to be independent of k 
and CT. 

Using the nonequilibrium Green's function technique, 
the spin-dependent current through the left ferromag- 
netic reservoir can be expressed as [48,49] 

ie f de 

lLa = T —Tr(ff,rL{G<(e)+fL[ &■■(£) -G^(e)]}), 



where Tr is the trace in the spin space; ct^ is a 4 x 4 matrix 
with Pauli matrix cr^ as its diagonal components; C'^'^is) 
are retarded, advanced, and lesser Green's functions in the 
generalized 4x4 Nambu notation. 

G\t,t) = -iOit - t ){{^it),^''(t)}). (7) 



where Hf is the Hamiltonian of the left and right ferro- 
magnetic electrodes 

Hf= J2 ^^v.kaC'^^ka^v.ka- (2) 
v=L,R ka 

Here, c^^ f^^{Cv,ka) is the creation (annihilation) operator 
in the lead v with energy £„ ,(;cr • represents two super- 
conducting reservoirs with chemical potential /Xg = 0 and 
the energy gap A, 

^s= XI ll''y-i<^''l,ka^y-'<''+H^^''l,kA,-ki+^-^-'^- 

Y=l,2 ka k 

(3) 

^DQD in Equation 1 denotes the DQD Hamiltonian 

^^DQD = X SidJ^dia + tc 'Y^{dl^d2a + h.C), (4) 

in which d\^ (dia, i = 1, 2) represents the creation (annihi- 
lation) operator of the electron with energy in the dot /; 
tc is the coupling strength taken as a real parameter The 
last term, H-x, in Equation 1 corresponds to the tunneling 
Hamiltonian between the DQD and four leads, 

= T.^T%clkJi.+h.c^+Y,{Ty^,kA''+h.c.). 

kiva kiycr 

(5) 

where the tunneling matrix elements between the DQD 
and two ferromagnetic leads are T^j. = ITliIc"^/*, T^j. = 
|7L2|e-''*/4, = |rRi|e-''^/4e-'*'Ki/4, and = 

|7R2|e''^/^e'*'R2/4'. The phase shift due to the total mag- 
netic flux threading into the AB ring is assumed to be 
(p = 2:/t(<I>l + <I'r)/0o with the flux quantum (po = h/e. 
The phase factor (p^i comes from the RSOI in dot i, which 
is tunable in the experiments [46,47]. T^'^^ = 7si(rs2) 
as the tunneling coupling between the DQD and two 



G^(t,t) = i{^Ut)^(t)), (8) 

with the vector ^I*^ = id\^, di^, d!^^, (^24, )• 

After some algebraic manipulations, the spin-depen- 
dent current can be derived from Equation 6: 

lLf = lf ds[ T^^{s)if^ -A) + -7g) 
+ rLK(e)(/L -/r) + Tf'isWi. -/s)] , 

(9) 

lLi = ^f ds[ T^^{s){fi, -Jo + rf^^£)(/R -7l) 
+ r[«(e)(7R -A) + (e)(/s -A)] , 

(10) 

in which T^^ and T^^^ are the spin-dependent AR and 
CAR coefficients, respectively. T]f' represents the single- 
particle tunneUng through FL-DQD-FR or FR-DQD-FL. 

OS 

Ta corresponds to the probability of the quasip article 
tunneling among two superconductors and the left ferro- 
magnetic lead./L A)i/r A)' ^''^^ /s '''^fi Fermi-Dirac dis- 
tribution functions. The derivation of the spin-dependent 
current is minutely given in the Appendix. 

Since we mainly focus on the AR process at zero temper- 
ature limit and set \eV]J\ = \eV^\ < A, will vanish. In 
the case of cVl = cVr, the current from the quasiparticle 
tunneling through FL-DQD-FR or FR-DQD-FL becomes 
zero; as a consequence, the AR dominates the transport 
through the four-terminal AB interferometer 

Results and discussions 

In the following numerical calculations, we mainly eluci- 
date the spin-dependent AR process in the four-terminal 
AB interferometer with the coherent indirect coupling 
and the RSOI. We take e = /j = /cb = 1, and set A = 1 
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as the energy unit. Throughout the paper, the symmet- 
ric couplings with TJ' = = Ts = T = 0.2A and 
I-PlI = I-PrI are considered as a typical case. 



Conductance 

Because we mostly investigate the AR within the super- 
conductor gap, in the limit of zero bias Vl = Vr — > 0, the 
spin-related AR and CAR conductances have the forms 



-AR 



2e' 



and 



gCAR^ 



2e^ 



rCAR 



(£f). 



(11) 



(12) 



It is well known that a DQD system with the maxi- 
mum coupling I a I = 1 has already been investigated. 
Indeed, such case is very special, and most experimental 
conditions correspond to la] < 1; as a result, a charac- 
terizing the coherent indirect coupling between two QDs 
via two ferromagnetic electrodes is introduced (see the 
Appendix), la] < 1 comes from the various factors, such 
as imperfections in the ferromagnetic reservoirs produc- 
ing the destructive quantum interference, the geometrical 
structure of the system, and so forth. 

Let us begin with the case of 0 = 0 and cp^ = jt/2; for 
the different coherent indirect coupling a, Figure 2 shows 
the total AR conductance (G^ = Ga^^^'' + Ga^^''^-' and 
G^^ = Ga^^^^^ + Ga^^^^^^) as a function of Fermi energy 
£f for parallel (P) and antiparallel (AP) configurations. In 
order to gain the clear physics, the Hamiltonian //dqd 
is diagonalized, and two energy eigenvalues are given as 
E± = + £2) ± - £2)^ + 4t2]; thus, when the 

Fermi level coincides with the £+ and £_, the resonant 
AR occurs and two peaks of AR conductances are located 
around the level E± as illustrated in Figure 2a, b, c. For 
a = 0, it is clearly seen that is always equal to G^ in the 
P arrangement (/"l = Pr = 0.4), and the magnitudes of 
two peaks are equal. However, for the case of the AP con- 
figuration {Pi^ = -Pr = 0.4), G^P # G^P appears when 
a = 0. Because the ferromagnetic leads have majority and 
minority electrons, the AR and the CAR are governed by 
the minority electrons for P configuration; thus, the AR 
and the CAR do not contribute the spin-polarized trans- 
port. For AP alignment, although the AR cannot produce 
the spin-polarized current, G^^ 7^ G^^ emerges due to 
the CAR process, in which the CAR is governed by the 
majority electrons. This leads to the appearance of the 
spin-polarized conductance. Since two dots are indirectly 
coupled via two ferromagnetic leads, which is reflected in 
the nondiagonal coupled matrix elements (see Equations 
15 and 16), a = 0 implies that the off-diagonal matrix 



S 0.04 - 



^ 0.06 




Figure 2 The AR conductance versus Fermi energy for P and AP 
configurations, (a) a = 0, [b) a = 0.5, and (c) a = 1 .0. Other 
parameters are £1 = £2 = 0, fc = 0.5, 0 = 0, and (^r = jr/2. 



elements vanish due to complete destructive quantum 
interference; thus, two dots are totally decoupled through 
two ferromagnetic leads. The AR (or the CAR ) can hap- 
pen only through QDl and QD2, respectively. This leads 
to the conductance G^ = G^ for the P arrangement and 

the equal height of two peaks (G^^ or G^^) for the AP 
configuration. 

We also notice that both G^ # GJ and Gf # Gf^ 
occur with the increase of a for P and AP configurations; 
thus, G^ — G^ is nonzero ata ^ 0, which means the occur- 
rence of the spin-polarized AR for P configuration in the 
presence of the RSOI and the nonzero parameter a. As a 
matter of fact, we have found that G^ = G^ is indepen- 
dent of the parameter a for P configuration in the absence 
of the RSOI, which is not shown here. In comparison with 
the case of a = 0, the symmetry of AR conductances 
with respect to the Fermi energy is significantly broken 
when a 7^: 0. It is noticeable that amplitudes of conduc- 
tance peaks near the level £+ decrease, and the magnitude 
of the right peaks is smaller than that of the left ones. 
In addition, the positions of peaks for AP alignment are 
also shifted with a = 1. These results indicate that the 
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coherent indirect coupling and the RSOI play an impor- 
tant role in determining the feature of the AR conductance 
spectra. 

To elucidate better the properties of the AR under P and 
AP configurations, in analogy with the conventional tun- 
neling magnetoresistance (TMR) effect of ferromagnetic 
tunnel junctions, the ARTMR is introduced and defined as 

[(Gf + Gf)-(G^^ + Gl)] 



ARTMR = 



(Gf ■ 



Gf) 



(13) 



In Figure 3, we present the 4> dependence of ARTMR for 
different a. The oscillation period of the ARTMR versus 
magnetic flux (pis lit, and the sign of the ARTMR does not 
change. It is interesting to note that the convex shape of 
the ARTMR at (/> = 2«7r (« is an integer) relies on the sign 
of the coherent indirect coupling parameter a. In compar- 
ison to the case of la] = 0.5, the magnitudes of ARTMR 
are considerably increased for |a| = 1.0. This is because 
the reduction of the destructive interference results in the 
enhancement of ARTMR. 

As we know, the RSOI can induce the spin precession 
and may even cause the inter-dot spin-flip effect. Accord- 
ing to [26,27], the spin-dependent phase factor due 
to the RSOI can be expressed as = <pRi — (p^2 = 
Pm* (Ll — L2) If?, where /J is the RSOI strength, m* is the 
electron effective mass, and Li is the length of dot /. ^r is 
tunable in experiments. It can reach 7r/2 easily or can be 
larger experimentally [27]. In order to explore further the 
influence of the coherent indirect coupling and the RSOI 
on the ARTMR, the ARTMR as a function of the parame- 
ter a for different ^r is shown in Figure 4. We can see from 
Figure 4 that ARTMR versus a exhibits the nonmonotonic 
features, and there exists the crossing point at a = 0. 
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Figure 3 ARTMR versus the magnetic flux ^ witli different a. 

Other parameters are £1 = £2 = 0. tc = 0.5, and <pR = :t/2. 
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Figure 4 ARTMR versus tlie coherent indirect coupling a with 
different ^pj^. Other parameters are £1 = £2 = 0, fc = 0.5, and 

0 = 0. 



Since a = 0 means that the coupling off^-diagonal terms in 
Equation 16 are totally suppressed, as a consequence, the 
AMTMR is independent of the RSOI (see the Appendix). 
When <pR is relatively small, this corresponds to the weak 
RSOI strength; thus, the variation of the ARTMR with a is 
not smart (solid line and dashed line). However, the evo- 
lution of the ARTMR is very remarkable as ^r increases 
(dotted line and dash-dotted line), while the ARTMR first 
increases and decreases with the increases of a, even the 
negative ARTMR also emerges, which corresponds to a 
spin valve effect in the AR process. This reflects that the 
strong RSOI gives rise to the significant variation of the 
ARTMR. We also observe that the maximum and mini- 
mum values appear in the curves of the ARTMR. These 
demonstrate that the optimal ARTMR can be tuned by 
means of the external parameters. 

Spin-dependent current 

Above, we analyze the properties of the AR conductances. 
In the following discussions, we will explore the spin- 
dependent current in the AR process with the help of 
the current formulas (Equations 9 and 10). To gain a 
full physical picture on the DQD levels' influence on the 
spin-related current. Figure 5 displays the images of the 
spin-polarized current /j = /l^ — /l; as a function of 
the energy levels ei and £2 of the DQD. The blue regions 
correspond to zero current, namely, /l^ = /l; in these 
regimes. In the diagram, it is found that the spin-polarized 
current is symmetrical about the line of ei = £2 and is 
asymmetrical with respect to the line of e\ = —£2 as illus- 
trated in Figure 5a, b. It is interesting to note that one 
level is aligned to the Fermi level, and the other is far from 
the Fermi one (off^resonance). Is is relative small. This is 
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Figure 5 Images of spin-polarized AR current current as a function of QD levels ei and e2 (Color on line), (a) a ■ 

Other parameters are 0 = 0, ijsp, = 0, fc = 0.5, and Pl = Pr = 0.4. 



-0.5 and (b) a = 0.5. 



because one QD is in the on-resonance state and the other 
is in the off-resonance state. When both ei and £2 are close 
to the Fermi level by tuning the gate voltage, the maxi- 
mal Is appears since DQD is in the on-resonance states. 
We also observe that, for a = 0.5 and a = —0.5, the 
spin-polarized current shows the opposite feature, which 
is a swap effect originating from the different sign of the 
parameter a. This indicates that the sign of the coherent 
indirect coupling parameter has a remarkable impact on 
the spin-polarized current. 

As we know, when ferromagnets are fully polarized, 
two ferromagnets become half metals where all elec- 
trons have the same spin. AR is usually suppressed at 
the ferromagnet/superconductor interface. However, AR 



still can occur, and the pure spin current can be gen- 
erated in the present system. For Pl — —Pr = 1-0 or 
Pl = —Pr = —1.0, i.e., two ferromagnetic leads become 
two half metals; the normal AR vanishes due to T^^ (see 
Equations 21 and 25). However, CAR dominates the trans- 
port through the four-terminal AB interferometer for AP 
alignment. As a consequence, we can obtain the pure 
spin current via CAR and two half-metal reservoirs. Thus, 
this device may be used as a pure spin-current injector 
even in the absence of the RSOI. In Figure 6, we also 
depict AB oscillations of the spin current for different 
a. For the case of a = 0.5, the magnitudes of the res- 
onant peaks and valleys are enhanced with the increase 
of the RSOI strength, and positions of peaks and valleys 




Figure 6 The spin current versus the magnetic flux 0 for different ^r. (a) a = 0.5 and (b) a = 1 .0. Other parameters are £1 = E2 : 
Pl = -Pr= 1.0. 



: 0, fr = 0.5, 
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are also shifted to the left, as illustrated in Figure 6a, b. 
Since the RSOI gives rise to an extra spin-related phase 
factor (pJ^ (see Equation 16), the curves of the spin cur- 
rent versus magnetic flux (p move towards the left with 
the emergence of the RSOI phase, and the shifted mag- 
nitude of peaks (or valleys) is equal to ^or as shown in 
Figure 6a, b. Physically, the increase of cp^ corresponds 
to the strong RSOI, which also favors the CAR process 
and the generation of the large spin current. When the 
DQD is fully coupled via two ferromagnetic reservoirs 
(a = 1.0), in comparison with the case of a = 0.5, it is 
noted from Figure 6b that not only the positions of peaks 
and valleys are altered, but also the amplitudes of those are 
remarkably enhanced. This originates from the fact that 
the reduction of the destructive interference enhances the 
spin current for the case of a = 1.0. These results indi- 
cate that the variation of the spin current is sensitive 
to the parameter a and the strength of the RSOI, and 
the interplay between them determines the nature of the 
spin current. 

Conclusions 

In this paper, we have analyzed the AR of a four-terminal 
AB interferometer containing a coupled DQD with with 
the RSOI and the coherent indirect coupling via two 
ferromagnetic leads. The formulas of the transmission 
coefficients are derived based on the framework of the 
nonequilibrium Green's function technique. For P con- 
figuration, the spin-polarized AR can occur, stemming 
from the RSOI and a nonzero coherent indirect coupling. 
On the contrary, for AP configuration, the spin-polarized 
AR always happens because of the CAR mechanism. 
Under the introduction of the ARTMR, we find that 
the sign of the ARTMR versus the magnetic flux keeps 
invariable for different parameter a, but the convex shape 
of the ARTMR depends distinctly on the sign of the 



parameter a. With the increase of the RSOI strength, 
the ARTMR versus the parameter a exhibits the more 
significant nonmonotonic features, and there exist the 
extreme values in the ARTMR plot, even the negative 
ARTMR also emerges. Since the energy levels of the DQD 
can be manipulated via the gate voltage, we can obtain 
the optimal spin-polarized current. A pure spin current 
can be generated via the CAR and two half-metal leads. 
Moreover, the strong RSOI and the reduction of the 
destructive interference {a = 1) favor the enhancement of 
the spin current. Thus, this device may become an effec- 
tive spin-current generator, and the pure spin current is 
tuned in terms of the magnetic flux, the RSOI strength, 
and so forth. These results offer the ways to manipu- 
late the spin-dependent transport via the four-terminal 
AB setup. 

Appendix 

In this appendix, we present the derivation of the current 
formulas in detail. 

Let ^"^(e) and G'^(e) denote the retarded Green's func- 
tion of the DQD without and with the coupling to the 
external reservoirs. In the Nambu space, ^"^(e) can be 
given as 

0 -tc 
s + si 0 

0 S - S2 

tr 0 



[/(£)]" 



0 



V 0 



0 

tc 

0 

£ + £2/ 



(14) 



Based on the following Dyson equation, the retarded 
Green's function of the system can be written as 
[G^(e)]-i =[/(£)-! - -£'], in which = E[ -h 
Eg + Egj + Eg2- The lesser Green's function G^(e) = 
G^(e)E<G^(e), where G^(e) =[G'' (£)]■•■ and E< = E< -h 
Eg -|- E^^ -|- E^. In the wide-band limit approximation, the 
retarded self-energy can be derived from the definition 



Ef = 



I 



0 
0 



■0/2 



J4>li 



If^ 2t 

0 

nL 

0 



0 



pL pL Q-i4>l1 



0 

pL 



(15) 



y — 



/ 



0 
0 



0 



/2 



0 

pR 
^2t 

0 



0 



2| 



(16) 
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— 



:rsi = --ripi(e) 



-'S2 



where PH^ = lit Y.^^ 7^'^* ^fuPva with py„ being the 
density of states of the spin a band in the lead v. We 
calculate the tunneling matrix element by means of the 

Bardeen's formula, i.e., T^'^ = j^[ifv,k(x)^ i^lji^* (x) 

ff*Cr)'^fv.kCr)'\ dS/2m^, where We is the effective mass, 
S is the region of the integration, iry^j^(f) is the wave func- 
tion of evanescent mode of the lead v, and ijf^ is the wave 
function of an electron localized in the QD /. Considering 

the propagation of electrons in the reservoir v, this prop- 
agation process (the wave number dependence of T^'^) 
induces the coherent indirect coupling via the reservoir v 
between two QDs, which is characterized with the param- 
eter . We assume that (X,, Y,, 0) is the center position of 
the ith. QD, Xi = X2 = Xd and L = \Yi - Y2\. thus, a" 
is given by a" = a ~ {IX^f l{L^ + (2X£,)2]J based on 
[21]. We find la] < 1 and decreases with L, and |q;| = 1 
corresponds to i = 0. We define a = ^ua/^/'^ua^iia' 
in which Tj'^^ = F^^, = V^^. With the definition 
of the spin polarization = (F^^ - FJ^)/{FJ^ + FJ^) in 
the lead v, the tunneling matrix element can be written as 

r!^ = r[(i +p,) and r;; = r;(i - p,) with r; = (fj^ + 

F^^)/2. r^=i,2 = 27T Zka \T%\^Ny„, Nya is the density 
of states when the superconductor is the normal state, 
and Py (s) is the modified BCS density of states py (e) = 

Vs-' — A-' i\l A'-—B^ 

(p^ = (p^i — (p^2> the spin-dependent phase factor is given 
hy (pa = 4> + 2ct(^r. We mainly take account of the case 
of the symmetric coupling between two superconducting 
electrodes and DQD, that is, Fy = F^. According to the 
fluctuation-dissipation theorem, the lesser self-energy can 
be given as E< = - EJ) and T,<^ = iFyFgyPy, 

where py = pis) = (lel/Ve^ - A2)6>(|e| - A).f„ and Fy 
are, respectively. 



/ 1 


_ A 

£ 


0 0\ 


_ A 

£ 


1 


0 0 


0 


0 


0 0 




0 


0 0/ 


/O 0 


0 


0 \ 


0 0 


0 


0 


0 0 


1 


_ A 

£ 


vO 0 


_ A 

£ 


1 / 



(17) 



(18) 



//„(£- eK) 0 0 0 \ 

0 Me + eV^) 0 0 
0 0 Ms-eV^) 0 

V 0 0 0 /v(e-l-eVv)/ 



//,(£) 0 0 0 \ 
0 fs(e) 0 0 
0 0 Us) 0 
V 0 0 0 fsie)/ 



(20) 



inwhich/i,(e-eV'v,) — fv,fv(e+eVv) = /„, and/s(£) are the 
Fermi distribution functions. By substituting Equations 
15 to 20) into Equation 6, we can obtain the spin-related 
current as shown in Equations 9 and 10. The AR (CAR) 
coefficients {T^^ and T^^^) and the probability of the 

quasiparticle tunneling {T^^ and T^^) can be calculated 

as 

+ r2\rL^|G^2(e)l' + r2\r2\|G^4(£)|2 



+ 2aT\^^T\^T\^Re{e-'*I^G\^{e)Gl^{e)] 



(21) 



+ 2«rL^yr^i?e[ e-^-I^G\^{s)Gl^{s)] 

+ 2«yrL^rL^rf^rR^i?e[e«*-*'''/2G^4(e)G^i(e)] 

+ 2ocFlJT\.T\.Re[ e'^I^Gl^{e)Gl^{e)\ 



+ 2«y rL rL rf rf , i?e[ e'^^+^''^Gl^{e)Gl, (e)] 



+ 2«rfy r2\i?e[ e-'*/' GJ2 (£) GI3 (e)] 

+ 2xxT\^^T\T\Re[ e-'^" I^G\^{s)Gl^{e)\ , (22) 

+ rf^rL^|G^i(e)p + rL^rf^|G^3(e)p 

+ 2aF\^/rf^^Re[e'l'''^^Gl3is)GUs)] 

+ 2ayr\^r^^rf^r^^Rel e'^1'''+^^'^Gl3(8)Gl,(e)] 



+ 2ay rL^rL^rf^rf^i?e[ e'^'t"-'^^/^G[3(e)Gl3(e)] 



+ 2aF^^^r\^r^^Re[e''^^^Gl^(s)Gl^{s)] 
+ 2aF^^^Ff^F^^Re[e'l"^^Gl3{8)Gls{s)] 
+ 2aFf^^/T\T%Re[ e'-O'^Gl, (e)G?i (e)] , 



(19) 



(23) 
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T^\e) =rsp{r^^[ |Gii(e)|^ + \G\:,(s)\' + \Gl^ie)\ 

+ |G^3(e)|2 + |G^4(e)|2] 

+ 2«yrL^i?e[e'*/2G^2(e)G2i(e)] 

+ 2a^r\^r^^Re[e'1>^''Gl3(e)Gl,(s)] 

- 2r\^ jRe[ Gl2(e)G?i(e) + G[^is)Gl^ie)] 

- 2r^^ -Re[ Gl2(s)Gl3(e) + G^3(e)G^3(e)] 



■:^fT-'L r\r^r /„m2 , |/-.r ,-^^12,1/^1 /^\|2 



- 2ayrJ;^rL^7?e[e''^/'G^i(e)G^i(e) 
+ e'*/2G54(e)G^i(£)] 



- 2a^r\^r^^Re[e''^^^G';i2(s)Gliie) 
+ e'*/2G^3(e)G^i(e)] }, 



(24) 



+ 2arL^7rL^i?e[e-'*/2G^3(e)G?2(e)] 



+ 2ayr\^ r_L, rL^7?e[e-'^G^3(£)G?2(e)] 
+ 2«V^it^U^2tr^^^e[ G^i(e)G|2(£)] 

+ 2arL^yr[^7?e[e-'*/2^^^(e)G?2(e)] 



(25) 



+ 2ar[^ yrf^ rf^;?e[ e'^-^^gr^ (e)G^2(e)] 

+ 2ayrL^rL^rf^rf^7?e[e'<^'-«/2G^3(e)Gl2(e)] 



+ 2ayrL^rL^rf^rf^i?e[e'(^"+«/2G^3(e)G?4(e)] 
+ 2arf^yrLl^i?e[e'^/2G53(e)G^4(e)] 
+ 2arf^yrLl^i?e[e'^/2Gr^(g)G?4(6)] 
+ 2arL^yrK^i?e[e'^''/26^3(e)6?4(e)] 

(26) 



+ 2«rf^yr}^i?e[e'*-/2G^2(£)G|4(£)] 



+ 2«Vri\rL,rf,r«,i?e[e'('^"-«/2G^4(£)G^4(e)] 



+ 2ar2\yrL^rL^i?e[e'*/2G^4(e)Gt4(£)] 
+ 2«rL|VrRrf^i?e[e-'*''/2G^(e)G^4(e)] 
+ 2«rL^yr^^7?e[e-'*-/2G^4(e)G^2(e)] > 



(27) 



rQ^(e) =rsp{rL^[ |G^i(e)|2 + |G^2(£)I' + \Gl3(s)\^ + IG^4(e)l'i 
+ rL^[ |G^i(e)|2 + |G^2(£)I^ + |G43(e: 
+ 2a^/Y\Y\Re[e-'1'/^Gl^(s)Gl2(e)] 
+ 2a^r\^r^^Re[ e-^^^Gl2(s)Gl2(e)] 
+ 2a,/rL 7?e[e-'^/2G^3(e)G^2(£)] 



i;^ 2; 



+ 2ayrL^rL^i?e[e-''^/2g^^(g^(5a^^g^j 

- 2r}-^ |i?e[ G^2(«)G?2(£) + G^4(e)G|2(£)] 

- 2rL^ jRe[ Gl2(s)Gli(s) + G^(e)G|4(£)] 



-2«yrL^rL^7?e[e-''^/2G^2(e)Gi2(e) 
+ e-'^^^Gl^{E)Gl2is)] 
- 2a^r\^r^^Re[e-'^^^Gl^{e)Gl2ie) 
+ e-'*/2G^3(e)G^2(e)]}- 



(28) 



Thus, we can investigate the quantum transport 
through our model system based on the above-mentioned 
equations. 
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